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Abstract 



Spin relaxation in a microscopic model of spin maser is studied 
theoretically. Seven qualitatively different regimes are found: free 
induction, collective induction, free relaxation, collective relaxation, 
weak superradiance, pure superradiance, and triggered superradi- 
ance. The initiation of relaxation can be originated either by an 
imposed initial coherence or by local spin fluctuations due to non- 
secular dipole interactions. The Nyquist noise of resonator does not 
influence processes in macroscopic samples. The relaxation regimes 
not initiated by an imposed coherence cannot be described by the 
standard Bloch equations. 
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1 Introduction 



Several experiments have been recently accomplished studying the peculiar- 
ities of spin relaxation in nonequilibrium systems of nuclear spins coupled 
with a resonator. Such systems include Al nuclear spins in ruby (AI2O3) 
[1], and proton spins in propanediol {C^H^Oq) [2-5], butanol (C^HgOH) , 
and ammonia (NH 3 ) [6]. Radiation processes in these systems occur at 
radiofrequencies, that is, in the maser frequency region. Similar processes, 
also occurring at radiofrequencies, exist in ensembles of electron spins (see 
discussion and references in [7,8]). The main attention in studying these 
processes has been payed to the possibility of coherent effects in spin relax- 
ation, the effects that resemble the well known optical superradiance [9-11]. 

For describing spin relaxation, one usually invokes the Bloch equations 
supplemented by the Kirchhoff equation for resonator. In addition, to solve 
this complicated system of equations, one resorts to adiabatic approxima- 
tion, when the resonator feedback field is proportional to transverse mag- 
netization. The latter approximation is, of course, unphysical for transient 
phenomena and could be used only at the final stage of relaxation. More- 
over, the Bloch equations as such, even being solved in a more realistic 
approximation [7], have the following generic limitations: (1) The role of 
dipole spin interactions is reduced to the inclusion of the spin-spin relax- 
ation time T2 ; (2) Coherent effects can be induced either by external 
magnetic field acting directly on spins or by imposing coherent initial con- 
ditions for transverse magnetization; (3) No self-organized coherence, that 
is called pure superradiance, can appear from a noncoherent state. These 
deficiencies do not permit to reach good agreement with experiment. 

To give a realistic picture of spin dynamics in nuclear magnets, one has 
to deal with a microscopic model. This kind of model has been treated by 
means of computer simulations [12-14]. However, great number of parame- 
ters in the model makes it practically impossible to analyse different relax- 
ation regimes by varying these parameters turn by turn. Here we present 
an analytical solution of the problem obtained with the help of a method 
developed in Refs. [15-17]. 



2 Relaxation Regimes 



The microscopic model for a system of nuclear spins is given [18] by the 
Hamiltonian 

1 N N 

H^^Hv-^BSi (1) 

i+3 »=1 
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with dipole spin interactions 

SiSj 3 ^Si^ij^J (^Sj^ij^j 



H -£ 

12 V r 3 

ij 



(2) 



where /i is a nuclear magneton and 

l 3 — I 1 3 I ' l 3 — 1 3i 



r ij 
r ij 



The total magnetic field 

B=H + H (3) 

consists of a constant external field Ho and an alternating field H of a 
resonator coil, 

i/ = # e z , H=He x . (4) 

The equations of motion, to be considered in what follows, will involve 
the average transverse spin component 

« = ^B^-iS?) (5) 
Iy i=i 

and the average longitudinal component 

s=^xm (6) 

iv i=i 

in which (...) implies statistical averaging. Since (5) is complex, we will 
need one more equation either for u* or for 

v = \u\. (7) 
For the resonator magnetic field we introduce the dimensionless notation 

in which 73 is the ringing width; uj , the resonator natural frequency; and 
Q is the resonator quality factor. The resonator is a coil of n turns of a 
cross-section area A . For the electromotive force Ef it is convenient to 
pass to the dimensionless quantity 



f=-2gf=f < X xwt. (9) 
nA h>y£ 
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For the functional variables (5)-(8), we can derive [16,17] the system of 
equations 

du 

— = i(cj Q - £ + i l2 ) u - i(-y 3 h + g)s, (10) 
ds % % 

-j t = 2 (7s/i + u* - - ( l3 h + C)u- 7 i(s - Soo), (11) 
dv 2 

_ = _2 72 t; 2 - i( 7 gfc + £)su* + i{ l3 h + C)su, (12) 

— + 2 73 /> + cu 2 / h(r)dr = -2a - (u* + «) + 73/, (13) 
at Jo dt 

in which 7l is a spin-lattice width; Soo , a stationary magnetization per 

spin, and 

Uo = — 72 = 1T' 

_ PV 2 ^ 72 
«0 = TT^T = T™? — > 

«73 73 

where p is the density of spins and r\ , a filling factor. The set 

f={£o,£,f} 

of stochastic variables models local spin fluctuations, £ being responsible 
for secular dipole interactions while £ and £* for nonsecular dipole inter- 
actions. The distribution of these stochastic variables is assumed to 
be Gaussian with a width 72 corresponding to a nonhomogeneous width. 
Equations (10)-(13) are completed by the initial conditions 

«(0) = «o, s{0) = z Q , v(0)=v o , h{0) = 0. (14) 

To solve the system of nonliner integro-differential equations (10)-(13), 
a method has been developed [15-17] combining the ideas of the multi- 
frequency averaging method, guiding center approach, and of generalized 
asymptotic expansion. Here we shall delineate the main steps of the method 
only in brief. 

2.1 Classification of functional variables 

Take into account that there are the following small parameters 

— <1, — <1, — <1, — <1. (15) 

72 ^0 ^0 ^ 

Also, the quasi-resonance case will be considered, when the detuning is 
small: 

|A| 

(A = lo-uj ). (16) 

^0 



As compared to the characteristic times 



T. ^ rp ^ rp* rp ^ 
1 = , J 2 = , J-2 = ~5 J 3 = , 

7i 72 72 73 

the spin oscillation time 

T = — « min^, T 2 ,T;,T s } (17) 

is small. 

If all small parameters in (10)-(13) tend to zero, then 

du , ds 

This shows that the variables u and h can be classified as fast and s 
and v 2 , as slow. 



2.2 Solution for fast variables 

The slow variables have the meaning of quasi-integrals of motion. Solve the 
equations (10) and (13) for the fast variables under fixed slow variables 

s — z, \u\ — V , 

treated as parameters in the corresponding solution 

u — u(z, v,t), h — h(z,v,t). 



2.3 Temporal and stochastic averaging 

The right-hand sides of the equations (11) and (12) are to be averaged over 
the small oscillation time (17) and over stochastic local fields £ . This is to 
be done after substituting into these right-hand sides the solutions for the 
fast variables found at the previous step. The resulting equations 

yield the solution 

z = z(t), v = v(t). 
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2.4 Definition of guiding centers 

The solutions z(t) and v(t) play the role of guiding centers for the slow 
variables. This role for the fast variables is played by 

u = u(z(t),v(t),t), 

h = h(z(t),v(t),t). 



2.5 Generalized asymptotic expansion 

Corrections to the guiding-center solutions can be obtained by expanding 
the functions 

u — u + . . . , h — h + . . . , 

s = z + . . . , \u\ 2 = v 2 + . . . 

about the guiding centers and substituting these generalized asymptotic 
expansions into the equations (10)-(13). In the latter, one has to separate 
out corrections of different orders which leads to a system of linear equations. 

This program has been accurately realized [16,17] yielding the following 
results. Considering the role of the Nyquist noise of resonator [19], it was 
found [20] that the relaxation parameter caused by the Nyquist noise is 
inversely proportional to the number of spins in the sample, 7^ ~ iV -1 . 
The value of this parameter is negligibly small for macroscopic samples with 
N ~ 10 23 . 

The solutions for the slow variables read 

Kt)| 2 = f IL) S ech 2 (tzh) + 2 p 2 ) 8 (t), (18) 

s(t) = ^nh( t -^-)--. (19) 
912 V r / g 

Here 

n- n ^(72 - 73) 2 , on x 

3 = °° Uj (T,-T»)» + A» <20) 

is the effective parameter of coupling between the spin system and the 
resonator. Approximately, 

gr ^ A(72 -7s) 2 



( 72 - 73 ) 2 + A 2 



Note that g = if 72 = 73 and A ^ . The effective radiation width 
7o , related to the radiation time r as 70 r = 1 , is given by the equation 



7o 2 = r 2 + (<n 2 )^ 



7 2 ~ 2 



(21) 
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in which t> = \u \, z$ = s(0) , and 

r = 72(1 + gzo)- 

The delay time is 

(22) 



7o 



70 + r 

Qualitatively different regimes of relaxation can be classified according 
to the values of two parameters, gz ans gvo . 

1. Free induction: 

g\z \ < 1, < gv < 1, 

to < 0, t « T 2 . (23) 

This is the usual case of free nuclear induction, with the maximal coherence 
imposed at t = . 

2. Collective induction: 

gz > -1, 5^0 > 1, 

to < 0, r < T 2 . (24) 

The coupling with the resonator is sufficiently strong to develop collective 
effects shortening the relaxation time r . 

3. Free relaxation: 

g\z \ < 1, v = 0, 

to < 0, r w T 2 . (25) 
This corresponds to an incoherent relaxation due to the local fields. 

4. Collective relaxation: 

gz > 1, v = 0, 

to < 0, r < T 2 . (26) 

Collective effects shorten the relaxation time, as compared to the previous 
case. 

5. Weak superradiance: 

-2 < gz < -1, v = 0, 
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to > 0, t w T 2 . (27) 
Here the delay time is positive; the process is weakly coherent. 

6. Pure superradiance: 

gz < -2, v = 0. 

t > 0, r < T 2 . (28) 

The coherence arises as a purely self-organized process started by local 
random fields and developed owing to the resonator feedback field. 

7. Triggered superradiance: 

gz < -1, gv > 1, 

to > 0, to < T 2 . (29) 

This is a collective but not purely self-organized process, since the relaxation 
is triggered by an imposed initial coherence. 

The term superradiance is used by analogy with the optical superra- 
diance [9-11], although, as is discussed in Refs. [17,20], there are many 
principal differences between the atomic and spin superradiance. 

The spin polarization, starting at z — s(0) , tends to 

s (t) ~ 2° (T 2 - r ) , t > t , (30) 

which shows that polarization reversal from negative zo to a positive s(t) 
occurs when r < T 2 . 

To illustrate the numerical values of the parameters occurring in the re- 
port, let us write down some characteristic quantities typical of experiments 
with proton spin systems [2-6], such as propanediol ( C 3 H 8 02 ), butanol 
( C^HgOH ), and ammonia ( NH 3 ). The spin of a proton is S = \ . 
The density of spins is p ~ 10 22 cm -3 — 10 23 cm -3 . The experiments are 
usually accomplished at low temperature T ~ O.lif in a magnetic field 
Hq ~ 10 4 G . The Zeeman frequency is u$ ~ 10 8 Hz , and the corresponding 
wavelength is A ~ 10 2 cm . The characteristic times are 

Ti ~ 10 5 sec, T 2 ~ 10- 5 sec, T 3 ~ 10~ 6 sec - lO^sec. 

The parameters in (15) are really small: 

^~10- 10 , ^~HT 3 , ^i~10- 3 , ^~10- 2 . 
7 2 luq lu uo 
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The relaxation parameter 77V due to the resonator Nyquist noise, as com- 
pared to the spin-spin relaxation parameter 72 , is 

™ ~ 10- 21 . 

72 

The picture obtained analytically is in agreement with the numerical 
simulations [12-14]. Some results of the latter are shown in Figs. 1-4, where 
K co h = Pcoh/Pinc is a coherence coefficient being the ratio of the coherent 
part, P co h , of the current power, P , to its incoherent part, Pi nc , and 
P z = —z(t) . These figures illustrate the qualitative changes of the processes 
under the variation of some parameters. For qualitative understanding, the 
absolute values of units in figures are not important, so they are not specified 
(for details see [12-14]). 

3 Conclusion 

The main results of the analytical solution for the problem of spin relaxation 
in nuclear magnets can be formulated as follows. 

1. Taking account of direct spin-dipole interactions is important for the 
correct description of relaxation processes. 

2. Nonsecular dipole interactions play the major role in starting spin 
relaxation, when no initial coherence is thrust upon the system. 

3. The Nyquist noise of a resonator does not influence the relaxation 
processes in macroscopic samples. 

4. There are seven qualitatively different regimes of relaxation, three 
of which are triggered by primary coherence imposed upon the system and 
four others are initiated by local dipole fields. 

5. The regimes that are purely self-organized, when no initial coherence 
is imposed, cannot be described by the standard Bloch equations. These 
regimes are: free relaxation, collective relaxation, weak superradiance, and 
pure superradiance. 
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Figure Captions 

Fig.l 

Coherence coefficient K co h , current power P , and the z -projection 
of spin polarization P z as functions of time, in arbitrary units, for two 
different coupling parameters, g\ and g 2 , with the relation g 1 / g 2 = 10 . 
The solid line is for g 1 ; the dashed, for g 2 . 

Fig.2 

The same as in Fig.l, for two different Zeeman frequencies, c^oi an d 

u 02 , related by the ratio u i/u 02 = 5 . The solid line is for u 01 ; the 
dashed, for CJ02 • 

Fig.3 

The same functions as in Fig.l, for different initial spin polarizations, 
with the relation 201/202 = 2 . The solid line is for z i ] the dashed, for 
Z02 ■ 

Fig.4 

The same functions as in Fig.l, for different initial transverse magneti- 
zations, with the relation foi/fo2 = 0.5 , when zo = . The solid line is 
for vqi ; the dashed, for v 02 ■ 
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